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Abstract 

Phonon polarization in a magnetic field is analyzed in 2D 
model. It is shown, that at presence of spin-phonon inter- 
action phonon possess elliptic polarization which causes the 
appearance of heat flux component perpendicular both to 
temperature gradient and magnetic field. 

PACS: 66.70. +f, 72.15.Gd, 72.20.Pa 

Recently a magnetic field dependence of heat conduc- 
tivity in dielectric crystal TbsGa^O^ has been experi- 
mentally established [T], [2|. It is so called Phonon Hall 
Effect (PHE) - a temperature gradient the has been 
measured in a direction, perpendicular both to a heat 
flow and magnetic field. This phenomenon is caused by 
spin-phonon interaction (SPI) of phonons and param- 
agnetic ions. The theory of this phenomenon was con- 
sidered in works [3], [4]. We can expect similar effect 
to be detected in dielectrics where molecules has ro- 
tary degrees of freedom. This effect is known for a 
long time in gases as Senftleben-Beenakker effect [5]. 
Generalization of this phenomenon on molecular crys- 
tals is considered in work [6]. The theory of Hall effect 
in ionic and molecular dielectrics for three-dimensional 
case is rather cumbersome. Therefore it is interest- 
ing to consider the theory in two-dimensional model. 
Probably, such theory can be applied to some experi- 
ments on quasi-two-dimensional crystals (films or sur- 
faces of three-dimensional substances). The spectrum 
and phonon polarization in presence of SPI was found. 
Expressions for these polar vectors demonstrate non- 
trivial symmetry, which is important for PHE theory 
construction. Transverse heat conductivity coefficients 
and nondiagonal density matrix calculations from |3J for 
two-dimensional case are briefly reproduced. In conclu- 
sion, we estimate the amount of effect and discuss the 
results. 

We consider dielectric consisting of light atoms and 
paramagnetic particles (atoms or molecules) with mag- 
netic moment M. In the two-dimensional case we as- 
sume that M is perpendicular to 2D-crystal plane. For 
simplicity we will accept that in a crystal cell there is 
one paramagnetic particle. The magnetic moment of 
rare-earth atoms is caused by electrons of /-subshell, 
the magnetic moment of molecules is proportional to 
their rotary moment. In both cases the spectrum of 
paramagnetic particles has a complex structure. But at 
low temperatures several bottom levels play the main 
role, and in simplest approximation we can replace the 



magnetic moment by pseudospin s„ (n is a cell num- 
ber). We will assume that T c < T < 6, where T c 
is a rotary degrees of freedom freezing temperature, 6 
- Debye temperature. In these conditions atomic os- 
cillations in crystal are caused by long-wave acoustic 
phonons. Thus all atoms in elementary cell oscillate 
with the same amplitude U„ and velocity V„, and the 
general motion of atoms in a cell creates the total or- 
bital moment [U„ x P„], r,n,e P„ = moV„, and mo is 
the total mass of atoms in a cell. Interaction between 
the moments of paramagnetic particles and the orbital 
moment of a cell we will describe by the following spin- 
phonon interaction Hamiltonian [7]: 



Hi = 5 ^(s„,[U„ x P„]). 



In this work we use a system of units where ks = 1, H = 
1. On the acoustic phonon wavelength scale there is a 
magnetization self-averaging, and we can replace oper- 
ator s„ with a value (s) = (s n ) ~ (M„) averaged over 
a crystal. We consider the interaction constant g as a 
phenomenological parameter. The value of g strongly 
depends on crystal field and it was estimated in numer- 
ous works [7]. 

Let us find the acoustic phonon renormalization 
caused by this interaction in a two-dimensional crystal. 
First of all, SPI changes a relation between velocity and 
an momentum of oscillating atoms (see [3]) 



dU n dH 



dt 



dP r 



= P„/m + #)xD„] (1) 



The oscillation equation in the momentum representa- 
tion is 

= D k u ks (2) 
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Here Ufe s is the normalized polarization vector, k is the 
phonon wave vector, s is a number of acoustic mode, 
and 

Df = Df - ie abc G c G c = 2ug (s c ) . (3) 

We see, that SPI adds an imaginary antisymmetric ten- 
sor to the symmetric dynamic matrix D^. In the zero 
approximation on SPI, the solution of ([2]) determines 
the dispersion law for two acoustic branches and two 
corresponding orthonormal polarization vectors. The 
eigenvalue equation is 



i 



To simplify, we assume crystal surface to be square- 
symmetric [8] 

Df = A^k 2 + A 2 S ab k 2 a + A 3 k a k b . (4) 

It's easy to verify, that in isotropic model (A2 — 0) the 
longitudinal mode has the energy u! k \\ — k\J A\ + A3, 
and the transverse mode energy is (uiq = k^/Ai) . In 
the general case we have 

uj 2 = X -{2A X +A 2 + A 3 )k 2 + s^R, s = ±1, (5) 



R 



(A 2 + A 3 ) 2 (k 2 x ~k 2 ) 2 + 4(A 3 k x k y ) 2 (6) 



For brevity, we omit the k dependence indication in 
evident cases here and below. In the general case, all 
of Ai are positive and nonzero. So, in real 2D-crystal 
the acoustic phonon spectrum is nondegenerated for all 
k in zero approximation, and lo 2 + > max(D™, Df), 
LU 2 _<min{D x k x ,Df). 

From © and ((H we see that R 2 and u)\ are in- 
variant with respect to inversion and following reflec- 
tions: ik x ,ky) > [k x , ky) , {k x ,ky 

) — > (~k x , k y ). For 
heat flow calculation we also need anisotropic velocity 
of sound expressions: 



duJs 1 S(OA 



A 2 + A 3 



+s-[{A 2 + A 3 ) 2 (k 2 a k 2 a+1 ) + 2A 2 k 2 a+1 ]}k a . 

If all of Ai has the same order of magnitude, then both 
velocities of sound has similar to k direction, and c ~ 
V^4 on the order of magnitude. 

Polarization vectors play a great role in the PHE 
problem. From dispersion equation ^ we have 



, ks - Df) C ks signk x , e 
and normalization 



v 



1) 



■nj\ 



ka - i-fe iC ks signk v , 



(7) 



c ks 2 = H s - Df) 



Df 2 = sR (cL - Df) > 



Eigenvectors (O and C ks are determined to phase, 
which depends on k and s. Let us take C ks — \C ks \. 
Then normalized polarization vectors are real and can 
be written down as follows: 



= s * signk 



y r > 



(8) 



--ka - Si 9 nk 

They also possess properties: 

(e+e_) = 0, (e+ x e_) z = C k+ C k -DfR, (e s fc) = 

.- n ,.J\k-\> »H.-DZ v ) 1 J\kv\> sU.-D'tf 
~ S * V ~P R h V ~P R ■ 

(9) 

From the last equation in ^ we can see that the up- 
per (lower) branch is approximately longitudinal (trans- 
verse). Components |(8|) changes their signs after inver- 
sion and corresponding components changes their signs 



after reflection from general axes as e ks are polar vec- 
tors. If we will turn k from <p — to 2ir smoothly, 
the vector e k - changes its direction by a jump when k 
crosses one of the axes. Vector efc+ oscillates near wave 
vector k direction, vectors e k + and are always per- 
pendicular to each other. 

Let us now discuss phonon renormalization in lin- 
ear approximation on SPI. We present the polariza- 
tion vector in linear approximation in the form of 
u s = ( e s + ^ e s)> an d rewrite dispersion equation <[2j) 
using this form. Then from the real part of obtained 
equation we can see that phonon spectrum and group 
velocity c ks — duj ks /dk are not renormalized (Slu s = 0). 
Imaginary part of dispersion equation determines the 
renormalization of polarization vector: 



\0J 



2^ab 



D ab ) 5e b 



-ie a b c G c e 



The most general form of the solution would have a 
form 5e b = iKe b s , , then we obtain 



K = 



([e k+ x B k -]G). (10) 



Thus, SPI leads to elliptic polarization of phonons, 
which is expressed in the imaginary addition to polar- 
ization vectors in zero approximation JH]). Using rela- 
tion (J3j) between G and g, for thermal phonons we can 
estimate the degree of elliptic polarization: 



K ~ g (s) IT. 



(11) 



As noted above, phonon spectrum and group velocity 
do not depend on SPI in linear approximation. But the 
form of heat flow density p in crystal can change in 
presence of SPI. However, it can be shown that, due 
to linear relation between velocity and momentum (Q]), 
the expression for heat flow in coordinate representation 
has exactly the same form as in [9] 
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ab 



d: 
13 « 



u?v b 



(12) 



but here V„ = P n /m + g[s x U„] (see. ([J)). Then 
we write lfl2|) in secondary quantization representation, 
average it over the system state, omit the anomalous 
averages (a ks a- ks >), (at ks ,a ks ), change symbols under 
the summation sign and obtain the following expression 
for the heat flow density: 



V1D 



ab\ 



(13) 

where p ss '(k) = (a ks aks') is phonon density matrix. In 
zero approximation (jT3j) has the usual form of energy 



flow density for phonon gas with n ks — (a ks ak s )- Sub- 
stituting polarization vector in linear approximation, 
we obtain (see [3]) 



(«y c ) = |E fc ( v .. / , 



^)K(k)w k ^u, k+ x 




(14) 
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If magnetic field and moment M are directed along 
z axis, temperature gradient VT is directed along x 
axis, and transverse heat flow is directed along y axis, 

then, substituting p |_ in the form of linear response 

p-+(k) = -iA x _ + (k) (VT) a into (JUJ), we have 




^ ± )K(k)u> k -Uj k+ x 



(15) 



where K(k) is determined in lflO|) . As K(k) ~ G, then 
for calculation lfl5|) in linear approximation, it would 
suffice to calculate A x 2 {k) in zero approximation. Ac- 
cording to [3] 



A x r= 
P9 



i^(k) 



F Ml (k)^(k)+F al ,(kK,(k) 

} 

AT fes (l + iV fes ) 



^ks Qpq 

2T 2 a 



(16) 
(17) 
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Here O p9 are the effective relaxation frequencies, p = 
ks, q — ks'. In general case O p9 , f2 gp , f2 pp have the 
same order of magnitude (but \Cl pq \ , \fl qp \ < fl pp ). In 
the two-dimensional case the longitudinal component 
of heat conductivity tensor is >c xx ~ T 2 f2 _1 , where c, f2 
are the average values of c p , ft pp . The substitution of eq. 
(JHl) and (fl7| into (fl5|) gives us the transverse compo- 
nent in the form of integral dependent on the relaxation 
frequencies ratio fl pq /fl pp . 

The product c x c y enter into the integral lfT5|). instead 
of (c x ) 2 in longitudinal heat conductivity, and it seems 
to disappear when averaging over k directions. How- 
ever, in ((TSj) there is K(k) ~ (ei x e 2 ) ~ D xy ~ k x k y . 



Therefore, under summation sign in (fT5|l there is an ex- 
pression invariant by inversion in reciprocal space, and 
averaging over k gives a nonzero result. We note, that 
the presence of ((ei x e%) G) points to determining role 
of the phonon elliptic polarization, which arises due to 
SPI. The integrand in l|15p becomes maximum when 
acoustic modes becomes at most close to each other. In 
2D crystal (when all A4 in |(4]) are the same order of 
magnitude) there are no such preferential directions. 

In contrast to 3D case [3] the PHE magnitude (flB")) is 
strongly definable 



77 y^X\l~^— + \l' : ^ L )K{k)uj k -UJ k+ x 



x(Cf 



fc+ 



CD 



i^(k)cJ(k) + F s , s (k)cJ,(k) 



UJ ks 



and can be calculated exactly. However, the problem 
definition has a model character and we will make the 
estimations only. We assume the collision frequencies 
fl pq and £l pp to be the same order of magnitude O. 
For simplicity we also accept a certain average value 
for velocity of sound c for both modes, and thermal 
phonons with frequencies about T play a main role in 
integral. In comparison to the longitudinal heat con- 
ductivity, the transverse heat conductivity has two ad- 
ditional multipliers: the degree of ellipticity ifTTj) and 



£l pq / [i-Oks — oJks')- They characterize the magnitude of 
nondiagonal elements of density matrix. The Hall angle 
has the following order of magnitude b 



3(f) 
T 



The last multiplier can be estimated from experimental 
value of heat conductivity coefficient. The first mul- 
tiplier depends on SPI magnitude and was discussed 
in works [5], [7], [TO] for ionic dielectrics. We suppose 
this ratio to have a comparable order of value in di- 
electrics containing molecules with rotational degrees 
of freedom. It can be even greater due to a higher de- 
gree of molecules anisotropy. Anyway, the transition 
from 3D to 2D case can make easier the observation of 
PHE. 
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